Search for the existence of circulating currents in high-T c superconductors 
using the polarized neutron scattering technique 
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We review experimental attempts using polarized neutron scattering technique to reveal the existence in high 
temperature superconductors of a long-range ordered state characterized by the spontaneous appearance of current 
loops. We draw a particular attention to our recent results (B. Fauque et al, Phys. Rev. Lett. 96, 197001 (2006)) 
that, up to now, can be explained only by the theory of circulating current proposed by CM. Varma. 



Since the discovery of high-T c copper oxide su- 
perconductors, a large variety of theoretical mod- 
els has been proposed to capture the intrinsic mi- 
croscopic nature of Cu02 planes, where supercon- 
ductivity develops. Several models postulate the 
existence of circulating currents (or current loops) 
in the Cu02 planes that could be responsible for 
the exotic electronic properties of these materi- 
als, in particular the existence of the mysterious 
"pseudo-gap" phase. In underdoped high-T c cop- 
per oxide superconductors, this phase is notable 
from its departure from the behavior of conven- 
tional metal jll2l.'{| H5l(il7| and it has been pro- 
posed that it could correspond to a long-range 
ordered phase associated with a new state of mat- 
ter, that could either co-exist or compete with 
superconductivity. 

In the staggered flux phase model [E] or d-wave 
charge density wave (DDW) model |9ll0lllll2j . 
CuC>2 planes are characterized by the appear- 
ance of staggered circulating currents flowing into 
square copper plaquettes (Fig.Q]a). In contrast to 
these models, where both translation invariance 
and time reversal symmetry are broken, the time 
reversal symmetry is broken in the circulating 
current (CC) phase proposed by CM. Varma, but 
translation invariance is preserved |13I14I15I16| . 
The order parameter of CC-phase is made of a 
set of staggered circulating currents embedded 
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Figure 1. Orbital currents in Cu02 plane (cop- 
pers: filled circles, oxygens: opened circles ): a) 
DDW state, b) CC-phase 0/, c-d) two domains 
of CC-phase 0//. The dashed lines indicate the 
size of the related Brillouin zone. 

in a given Cu square plaquette. Each current 
flows along a triangular path involving a copper 
site and two of its nearest neighbor oxygen sites 
(Fig.nb-d). Two distinct currents patterns have 
been proposed: the 0/-phase made of 4 triangu- 
lar loops (Fig. ^b) or the ©//-phase made of 2 
triangular loops (Fig. forming a "butterfly" 
centered on Cu site (Fig. QJc-d) 

Whatever its origin, the array of circulating 
currents in the Cu02 planes generates a magnetic 
field distribution: for the staggered flux phase or 
the DDW phase, one may speak about antiferro- 
magnetic (AF) orbital magnetism (with a prop- 
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agation wave vector at q=(7T, 7r) and q=0 AF- 
orbital magnetism in the case of CC-phases. Neu- 
tron scattering technique can be used to probe 
the existence of circulating currents in high-T c 
copper oxide superconductors. Neutron spin does 
indeed interact with any internal magnetic field. 
The magnetic intensity for scattering neutron is 
written |17|: 

(^)^ f K\<f\o.B(< l )\ l >f (1) 

where a is the neutron spin and \i > and \f > are 
the initial and final states of the neutron. B(q) 
describes the Fourier transform of the magnetic 
field distribution at the wave vector q. When the 
internal magnetic field originates from the cir- 
culating currents |8I12| . B(q) can be expressed 
in terms of the Fourier transform of the current 
density j(q), i.e. B(q) oc ^ c ^ ) 2 ><q . The neutron 
scattering cross section for polarized neutron then 
reads [T7|: 

(^)W oc \ <f \ a —p-\ t> \ (2) 

The neutron scattering intensity in the spin flip 
(SF) channel, /, has to fulfill a polarization selec- 
tion rule equivalent of the one for spin moments 
which in absence of chirality can be written: 

Ip//q. = + Ip//§± ( 3 ) 

where the unitary vectors q and q^ are respec- 
tively parallel and perpendicular to wave vector 
q in the scattering plane, and z is perpendicular 
the scattering plane (Fig. Eld). For unpolarized 
neutron and taking into account the charge con- 
servation constraint, the cross-section takes the 
simple form: 

The description of the current density could be 
quite complicated due to the nonzero size of the 
Wannier function, the mixing of different orbital, 
and many body effects when the current becomes 
substantial |8I12| . For sake of simplicity, it can 
be useful to associate an effective orbital mag- 
netic moment to a given current loop. The mo- 
ment is located at the center of the current loop 



and is aligned perpendicular to the CuC>2 planes, 
owing to the planar confinement of the currents. 
Therefore, in contrast to a spin moment on Cu 
site, the effective orbital moment should be much 
more spread in real space. As a consequence, one 
may expect circulating currents to give rise to a 
neutron scattered magnetic intensity decreasing 
faster than the squared Cu magnetic form factor, 
Fq u {<j), at large wave vector in agreement with 
Eq.El 

Let us consider, at first, the case of AF orbital 
magnetism. There were several observations of 
an AF ordering in the underdoped regime of su- 
perconducting YBa 2 Cu 3 6+x |f8H9l20| . When 
observed, this phase usually develops well above 
T c (close to room temperature) and displays an 
ordered moment of about 0.05-0.02 /ib- The po- 
larization neutron analysis reveals that the AF 
phase is dominated by moments aligned in the 
Cu02 plane, as in the insulating AF parent com- 
pound. Considering, on the one hand, that other 
high quality single crystals with a similar doping 
level do not show a similar order |2f| and, on the 
other hand, that impurity substitution out of the 
Cu02 planes (in the CuO chains) could induce an 
AF order at ~ 300K in optimally doped samples 
(T c =93 K) [22, it is likely that the observed AF 
order is not a generic property of the underdoped 
state and may be either remnant of the AF insu- 
lating state or induced by impurities or defects in 
the CuO chains. 

Beyond this AF spin order, Mook et al. |20I23| 
reported polarized neutron experiments suggest- 
ing the existence of a weak AF quasi-2D order, 
potentially characterized by: (i) magnetic mo- 
ments perpendicular to the CuC>2 plane, (ii ) a 
decrease of the neutron scattering intensity at 
large wave vector larger that F^ u (q). These ex- 
periments can be viewed as providing evidence in 
favor of the DDW model. However, it is worth 
pointing out that the estimated magnitude of the 
ordered moment is ~ 0.0025//B, i-e close to the 
experimental threshold of detection |22j . Further- 
more, the small number of studied Bragg reflec- 
tions was not sufficient to get solid conclusions 
concerning the evolution of the magnetic signal 
at large wave vector. Finally, the lack of repro- 
ducibility of the neutron scattering data severely 
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questions their relevance for cuprates in general. 
In addition, the study of the charge excitation 
spectrum by angle resolved photo-emission do not 
show any indication of a doubling of the unit cell, 
resulting from breaking of the translation sym- 
metry expected for the DDW phase. Therefore, 
there are no compelling evidences of circulating 
currents breaking translation invariance and lead- 
ing to an orbital AF state. 

Concerning the two possibles CC phases pre- 
serving translation symmetry proposed by CM. 
Varma |13ll4ll5flH] . the 6/-phase was not de- 
tected by earlier polarized elastic neutron scat- 
tering experiments |24I25| . The 8j/-phase has 
been later on proposed |15j to account for an 
angle resolved photoemission experiment |26| . 
Using circularly polarized photons, A. Kamin- 
ski et al |2fi| reported a dichroic signal in the 
Bi2Sr2CaCu208+<5 system indicating a time re- 
versal breaking symmetry in the pseudo-gap 
state. 

Recently, this 8//-phase has been searched 
with polarized neutron scattering using uniax- 
ial polarization analysis In contrast to the 
earlier inconclusive attempts, experimental evi- 
dences in favor of the existence of the CC 8//- 
phase were found |27|: a magnetic order, hidden 
in the pseudo-gap regime of high-T c superconduc- 
tors, has been observed. The study performed 
by Fauque et al. [23 covers a large part of the 
phase diagram of high-T c cuprates and show a 
remarkable reproducibility of the neutron results, 
that was missing in the previous studies devoted 
to DDW phase. Five different samples of differ- 
ent origin were indeed studied: 4 YBa 2 Cu 3 06+x 
samples in the underdoped (UD :T C < T C:max ) 
regime and one Yo.ssCao.ist^CusOg+x in the 
overdoped regime (T c > T c inax ). Hereafter, the 
samples are labeled as: UD54 (x=0.55, T c =54 
K), UD61 (x=0.6 - T c =61 K), UD64 (x=0.6 - 
T c =64 K), UD68 (x=0.75 - T c =68 K), OD75 
(x=l - T c =75 K). Polarized neutron data have 
been obtained on the triple axis spectrometer 4F1 
at the Orphee reactor (Saclay) (see Ref. [2Zj for 
the description of the experimental set up). 

In contrast to the DDW phase or the flux 
phase, the CC-phase does not break transla- 
tion invariance and the magnetic scattering in- 
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Figure 2. Temperature dependencies of the raw 
SF (o) and NSF (+) neutron intensity measured 
at q= (0,1,1) in sample UD64 (from [27|). 

tensity appears at the same wave vector as the 
nuclear scattering intensity. For the CC-phase 
8//, both intensities are superimposed since in- 
terference terms with the nuclear structure fac- 
tor cancel out due the existence of magnetic do- 
mains. The use of polarized neutron is then es- 
sential to identify the magnetic scattering gener- 
ated by the circulating currents. Most of the data 
of Ref. were obtained in a scattering plane 
(010)/(001)(Fig. EJ. In order to evidence small 
magnetic intensity, measurements have to be per- 
formed on the weakest nuclear Bragg peaks where 
the magnetic scattering is expected: the Bragg 
peak q=(0,l,l) offers the best compromise. So 
far, only the polarized neutron data in the nor- 
mal state are reliable. The effects specific to the 
superconducting state such as the evolution of the 
magnetic response associated with the circulating 
current through T c remains to be settled care- 
fully. 

Figure [21 a shows the raw neutron intensity 
measured at q=(0,l,l) for SF and NSF channel 
with P//q (sample UD64), for which the mag- 
netic scattering is entirely spin-flip and given by 
Eq. When cooling down from room tempera- 
ture, the NSF and SF intensities display the same 
evolution within error bars, until T<T maff , where 
the SF intensity increases noticeably, whereas the 
SF intensity remains flat. This behavior signals 
the presence of a spontaneous magnetic order be- 
low T mag . The normalized magnetic intensity 
Imag (see Ref. |27j) is ~ 1 mbarn (Fig[3Ja), i.e. ~ 
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Figure 3. a) Temperature dependencies of the 
normalized magnetic intensity (see text), mea- 
sured in the SF channel for 5 different samples: 
at q = (0,1,1) (o) and at q = (0,0,2) (•). b) 
Phase diagram as a function of hole doping |27| . 

10 -4 of the strongest Bragg peaks: it is therefore 
impossible to detect it with unpolarized neutron 
diffraction. 

The unusual magnetic order in a tempera- 
ture and doping range that cover the range 
where the pseudo-gap state is observed in 
(Y, Ca)Ba2Cu3C>6+x- For the 4 underdoped sam- 
ples, the full normalized magnetic intensity sys- 
tematically increases below a certain temperature 
T ma g whereas no magnetic signal is observed in 
the overdoped sample (OD75) (Fig. |3J|. The de- 
duced T mag , defined as the change of slope in the 
normalized intensity I mag , decreases with increas- 
ing doping and matches quite well the pseudo- 
gap temperature, T*, of the resistivity data in 
YBa 2 Cu306+ x |6I27| . Likewise, T mag decreases 
almost linearly and extrapolates to K close to 
the hole doping nh — 0.19, the end point of the 
pseudo-gap phase @] (Fig. |3Jb). 

To proof that the signal reported here is mag- 
netic, it should follow the polarization selection 
rule given by EqO The polarization analysis, 
shown in FigQ]a-c, indicates that the observed 
signal follows the expected selection rule. That 
unambiguously demonstrates the magnetic origin 
of the low temperature signal. Moreover, the ob- 



served magnetic peak at q=(0,l,l) is resolution 
limited along the (001) direction, showing that 
the magnetic order is characterized by long range 
3D correlations (at least at T=75 K)|27] . The 
polarized neutron data of Ref. [57] therefore pro- 
vide compelling evidences of the existence of a 
magnetic order hidden in the pseudo-gap state of 
high-T c cuprate superconductors. But can this 
magnetic order be ascribed to the 0//-phase ? 

For comparison of the neutron scattering data 
with what one could expect from the 0//-phase, 
one can express the magnetic intensity either in 
term of orbital magnetic moments (Ref. [27J) or 
in term of the Fourier transform of the current 
density: both approaches being equivalent at a 
qualitative level. Here, we propose to consider 
the second approach. Since the G/j-phase can be 
viewed as an array of butterfly triangular circu- 
lating currents centered on each Cu site (FigQJc), 
the Fourier transform of the current distribution 
takes the general form: 
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with Bi the unitary vector in Cartesian coor- 
dinates and G vectors of the Bravais lattice. 
According to this expression, the neutron in- 
tensity could be different from zero only for 
wave vectors along q={(0,K,L),(H,0,L),(H,H,L)}. 
Furthermore, since the current loops are stag- 
gered in each Cu plaquette, there is no uniform 
magnetic field in the Cu02 planes, this implies 
that the neutron scattering intensity vanishes for 
q=(0,0,L). The magnetic signal is indeed exper- 
imentally observed at q=(0,l,L) (Fig. 121 111 and 
vanishes at q=(0,0,2) (Fig. 02a) ). 

When the circulating currents are modeled 
with a set of thin wires along which flow a cur- 
rent with a given intensity, one obtains for the 
Aij terms in Eq. Ay(q) oc sin(ir(iH + 
jK)/2)/(i.H + jK). The terms Ay become ei- 
ther 1 or at a Bragg position. The term (3(q) 
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is equal to 1 for a monolayer system. Since 
YBa 2 Cu306+ x is a bilayer systems, current loops 
in the 2 CuC>2 planes could be either in phase 
(ip — 0) or out-of-phase (ip = tt/2), so that 
(3(q) oc cos(irLd/c + ip) (d=3.3 A stands for the 
distance between the planes within the bilayer.) 
Finally, for the scattering plane (010)/(001), the 
neutron cross-sections, given by Eqs. 12151 can be 
written for the circulating current model shown 
Figs-Hc-d as 



h 



'P//q 



z P//q 

P//<L. 



2 {TrLd/c + 



ip//2 







(6) 

(7) 
(8) 



The L dependence of the magnetic Bragg in- 
tensity measured in the SF channel for P//q and 
q=(0,l,L) is in agreement with Eq. |S|with ip = 
(Fig. EJe). One may speculate that this in-phase 
configuration could favor the tunneling of charge 
carriers along the c axis. However, in Fig0]a-c, 
the magnetic signal for Ip//q and for Ip//z do 
not simply match and a magnetic signal occurs 
also for P//qj_, in contrast to what is expected 
in Eq. 17181 This discrepancy could be related to 
the fact that the CC picture of Figs.^c-d assumes 
perfect Cu02 planes and, for instance, does not 
incorporate the dimpling of CuC"2 planes in the 
YBa 2 Cu 3 0g +x system. Further, spin degrees of 
freedom might also play a role in producing in- 
plane magnetic moments (necessary to account 
for a non-zero signal for P//qj_), by spin-orbit 
scattering as proposed in Ref. 2S for the CC- 
phase. 

Considering the reported neutron data , the 
observed magnetic order is qualitatively consis- 
tent with the existence of a circulating current 
state proposed by CM. Varma to account for the 
pseudo-gap |13ll4ll5lTfl] . The neutron data nev- 
ertheless suggest that a more realistic description 
of the Cu02 should be included, in addition to 
spin degrees of freedom. It is also important to 
emphasize that, if, on the one hand, the neutron 
data support the CC-phase model, on the other 
hand, alternative interpretations cannot be ruled 
out considering the limited amount of data avail- 
able. 
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Figure 4. Temperature dependencies of the nor- 
malized magnetic intensity, measured in the SF 
channel at q=(0,l,l) in sample UD61: a) for 
P//q, b) for P//z, c) for P//q_L d) Sketch of the 
scattering plane showing the three polarization 
directions discussed here, e) Sample UD54: In- 
dependence of the normalized magnetic intensity 
measured at T=70 K and q=(0,l,l) as compared 
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At that stage, one could consider, as an exam- 
ple, a ferromagnetic order with magnetic Cu mo- 
ment aligned perpendicular to the Cu02 planes. 
The magnetic intensity at q=(002) should vanish 
as well, whereas a signal should be observed at the 
Bragg reflection q=(01L), as observed experimen- 
tally. Furthermore, the L-dependence reported 
in Fig. ^e is dominated by the cosine square 
term, corresponding to " ferromagnetically" cou- 
pled planes. Since the studied Bragg reflections 
are at rather small |g|, one cannot distinguish be- 
tween an evolution at larger vector controlled by 
a term ~ for circulating currents or by a term 

~ (1 — ^r)F^, u (q) for ferromagnetic spins (for in- 
stance). However, this cannot account for the 
magnetic intensity observed for P//qj_ (Fig0]a- 
c). Another spin model could be a magnetic order 
involving staggered spins on oxygen sites which 
could also account for the observed magnetic scat- 
tering (see Ref. |27|). Thus, a spin order can- 
not be ruled out on pure experimental grounds. 
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The alternative spin orders have nevertheless two 
main drawbacks: (i) they should have typically 
been detected by local probe measurements such 
as NMR or [iSR, (ii) no theoretical model has 
been developed, so far, suggesting this could be 
relevant for cuprates. 

It is clear that more work is needed to de- 
termine the exact magnetic structure factor and 
the evolution of the magnetic intensity at larger 
wave vector. Further experiments will require a 
study of the magnetic signal with different ex- 
perimental set-up: (i) by using also a polarizing 
mirror to analyze the final polarization (ii) by 
performing the 3D polarization analysis 29 (us- 
ing CRYOPAD for instance) . Another advantage 
of CRYOPAD is the suppression of the magnetic 
guide field at the sample position which allows 
to study the magnetic response in the supercon- 
ducting state. Further, an external magnetic field 
could be applied since it should couple to spins, 
but should have a minor effect (if any) on the 
circulating currents. 

As a conclusion, we have seen that the neu- 
tron polarized scattering is a unique technique 
that has been used to search for the existence 
of circulating currents in cuprates yielding tiny 
magnetic field distribution to which neutron spin 
can couple. While the search for flux phase and 
DDW phase remains unsuccessful, there are now 
experimental evidences of a magnetic order hid- 
den in the pseudo-gap state of high-T c supercon- 
ducting cuprates, supporting the circulating cur- 
rent state proposed by CM. Varma in his theory 
of the pseudo-gap phase. 
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